Accepted to Icarus 



The Rotation of Janus and Epimetheus 

Matthew S. Tiscareno,"*^ Peter C. Thomas,^ and Joseph A. Burns"*^'^ 

^Department of Astronomy, '^Center for Radiophysics and Space Research, and 
^Department of Theoretical and Applied Mechanics, Cornell University, Rhaca, NY 14853 

ABSTRACT 

Epimetheus, a small moon of Saturn, has a rotational libration (an oscillation about 
synchronous rotation) of 5.9°ibl.2°, placing Epimetheus in the company of Earth's Moon 
and Mars' Phobos as the only natural satellites for which forced rotational libration has 
been detected. The forced libration is caused by the satellite's slightly eccentric orbit 
and non-spherical shape. 

Detection of a moon's forced libration allows us to probe its interior by comparing 
the measured amplitude to that predicted by a shape model assuming constant den- 
sity. A discrepancy between the two would indicate internal density asymmetries. For 
Epimetheus, the uncertainties in the shape model are large enough to account for the 
measured libration amplitude. For Janus, on the other hand, although we cannot rule 
out synchronous rotation, a permanent offset of several degrees between Janus' mini- 
mum moment of inertia (long axis) and the equilibrium sub-Saturn point may indicate 
that Janus does have modest internal density asymmetries. 

The rotation states of Janus and Epimetheus experience a perturbation every four 
years, as the two moons "swap" orbits. The sudden change in the orbital periods 
produces a free libration about synchronous rotation that is subsequently damped by 
internal friction. We calculate that this free libration is small in amplitude (< 0.1°) and 
decays quickly (a few weeks, at most), and is thus below the current limits for detection 
using Cassini images. 

Subject headings: satellites, dynamics — satellites, shapes — Saturn, satellites 
Running header: Rotation of Janus and Epimetheus 



1. Introduction 

Janus and Epimetheus, two small satellites of Saturn, share a remarkable orbital configuration. 
They move along the same orbit, but in contrast with other, much larger, bodies that are known 
to have small "Trojans" orbiting at their Lagrange points (including Jupiter, Mars, Neptune, 
Tethys, and Dione), neither Janus nor Epimetheus is overwhelmingly large compared to the other 



Janus is 3.6 times as massive as Epimetheus (Jacobson et al. 20081. Consequently, both Janus 
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Fig. 1. — The orbits of Janus and Epimetheus (with epicycles removed) in a frame rotating with 
their mass-averaged mean motion. The moons' radial excursions are exaggerated by a factor of 500, 
and the moons' sizes are exaggerated by a factor of 50, but otherwise the figure is to scale. Numbered 
points mark the moons' positions on 1) 2004 July 1, 2) 2005 May 21, and 3) 2006 September 9. 
Arrow lengths indicate motion accomplished in 100 days. 



and Epimetheus (to which we refer collectively as the "co-orbital moons") execute "horseshoe" 
trajectories about their mass-averaged mean orbital state (Dermott and Murray 1981a|b Yoder 



et al. 1983 Peale 1986: Tiscareno et al. 2006). Every 4.00 yr they reach their mutual closest 
approach and "swap" orbits (Fig. [T| — one moon's orbital rate (mean motion) slows down, which 
is to say that its mean distance from the planet (semimajor axis) increases, while the other moon 
does the opposite. The effect on Epimetheus is greater than the effect on Janus, in proportion to 
their masses. The most recent orbit swap occurred on 2006 January 21, at which time Janus became 
the inner satellite and Epimetheus the outer. The next orbit swap will occur on 2010 January 21. 

The periodic change in the orbital rate provides a unique opportunity to observe the rotational 



response of the moons. Tiscareno et al. (20061 made similar use of nature's obliging variation of 



a single parameter in a complex problem, observing morphological changes in the spiral density 
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waves raised in Saturn's rings by Janus and Epimetheus and connecting them with the periodic 
changes in the orbital rate. 

If one assumes that the co-orbital moon is rotating synchronously (keeping the same face 
always towards Saturn) as it goes into the orbit swap, the change in the orbital rate means that the 
rotation rate is no longer synchronized. Because Janus and Epimetheus are both significantly non- 
spherical, non-synchronous rotation does not continue indefinitely but results in simple harmonic 
motion (libration) about the most stable state in which the moon's long axis is pointed towards 
Saturn. This kind of motion, which results from a one-time impulse (the orbit swap) and then 
decays with time due to friction in the moon's interior, is a free libration. In principle, observing 
the free libration's rate of decay can yield insight into the moon's interior properties. 

There are, however, two other forms of libration that must be considered, both of which arise 
due to eccentricity of the moon's orbit. Firstly, because an eccentric moon moves faster at periapse 
and slower at apoapse, even a synchronously-rotating moon oscillates relative to the line connecting 
the moon and the planet (the optical libration). Secondly, the optical libration allows the planet 
to exert a torque on a non-spherical moon, resulting in a forced libration. The optical and forced 
librations are both in phase with the moon's orbit, with the moon's long axis pointing towards the 
planet at periapse and apoapse. To date, the only known moons for which a forced libration has 



been measured are Earth's Moon (Koziel 1967 Williams et al. 19731 and Mars' Phobos (Burns 



19721 [Duxbury and Callahan|[T98T| [T989t |Duxbury|[T989l [T99T| [Simonelli et~aL][T993| [Willner et al 



20081. With this paper, we add Epimetheus to that number. 



In the next section, we will begin by considering tidal de-spinning. We will then discuss the 
rotational free libration in the undamped (Section [s]) and damped (Section |4]) cases, and then 
the rotational forced libration (Section |5]) . In Section [6] we will compare our models with high- 
resolution imaging by the Cassini Imaging Science Subsystem (ISS). Section [7] presents a summary 
and conclusions. 



2. The Naive Case 

In its simplest form, the problem is characterized by the orbital longitude 6 and the libration 
angle ip (Fig. [2]). The orientation of the moon in inertial space \s 6 + ij). In the models presented 
here, we will make the simplifying assumption that the orbit swap takes place instantaneously, 
with the mean motion changing from ni to n2 at time t = 0. Since the moon was assumed to be 
rotating in a steady-state before the swap, ^{t < 0) = 0, the pre-existing rotation rate of the moon 
in inertial space is ni. The orbit swap then results in an impulse ipi = ni — n2. In fact, the orbit 
swap is not instantaneous, but takes place over a period of several months, making our estimate of 
'ipi an upper limit. 

The standard treatment of tidal de-spinning of a fast rotator is to consider a torque exerted 
by the central planet on a tidal bulge raised by the planet but offset from the moon-planet line due 
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Fig. 2. — The orientation ip of the moon's long axis is measured from the direction towards Saturn, 
while is the orbital longitude. All angles are positive in the counter-clockwise direction. Only in 
Section [5] is this picture complicated by non-circular orbits. 



to tidal dissipation, resulting in the expression (Peale 19991 



45 pn^R"^ 
76 fiQ 



(1) 



where ip is negative iff ^ > 0. The density p, mean orbital rate n, and mean radius R for Janus 
and Epimetheus are given in Table [T} 

The "quality factor" Q corresponds to the number of cycles it takes for rotational energy to 
substantially dissipate. A common assumption for a monolithic block of ice is Q ^ 100, though Q is 
likely much less for a fractured "rubble-pile" in which damping is made more efficient by increased 
internal friction. We will use Q ~ 10 to represent the latter case. The rigidity p is commonly taken 
to be ~ 4 X 10^ N m~^ for solid ice, but also decreases by an order of magnitude or more for a 
fractured body ( Goldreich and Sari]|2009 1 . 

The resulting values of ■0 are small — on the order of 1" dy^^ for solid ice, and some 100 times 
greater for a fractured rubble-pile. The de-spinning time for a fast rotator (say, Tpi = 100° dy~^) 
would thus be ipi/ip ~ 1000 yr (solid) or 10 yr (fractured). While these de-spinning times are 



Table 1: Parameters for tidal de-spinning 

R (km) m (10^** kg) p (g cm"^) n (°/dy) 
Janus 89^5 1.897 0.632 518.292 

Epimetheus 58.1 0.526 0.641 518.292 

Note: R taken from Table ksl below, while m and n are taken from|jacobson et al.| ||2008|l; p = 3m/4:nR^. 
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appropriate for the original damping of any primordial rotation of these moons, they are over- 
estimated for the case of the relatively small non-synchronous rotation induced by the orbit swap, 
ipi ~ 0.1° dy~^ (see Table [ij and recall from the beginning of this Section that this is an upper 
limit), which leads instead to ipi/i^ 1 yr (solid) or 4 dy (fractured). 

We would then appear to be in a felicitous regime for which solid and fractured bodies might 
easily be distinguished, as non-synchronous rotation induced by an orbit swap is quickly damped 
for the latter but may remain for the entire 4-yr horseshoe libration for the former. However, the 
above analysis suffers from two faulty assumptions: the first is that tidal de-spinning is controlled 
by the torque on the Saturn-raised tidal bulge, and the second is that de-spinning proceeds by 
slowing down a monotonic non-synchronous rotation. Because Janus and Epimetheus in fact are 
not spherical but significantly triaxial (Table [2|, the rotational dynamics are controlled rather 
by the torque on the moon's figure. Furthermore, because the moon's intrinsic figure does not 
change as the moon rotates, and because ■0i is not large, the moon executes a libration about 
the synchronous rotation state, with tidal effects appearing only as the mechanism by which the 
libration is damped. 



3. Free Libration: The Undamped Case 

We will now consider the limit of small oscillations, which our conclusions will show to be 
justified. In the undamped case, for small libration, the motion arising from the planet's torque 
on the moon's figure (modeled as a triaxial ellipsoid characterized by principal moments of inertia 
A<B <C) is described by (e.g., |Danby|[T988| Eq. 14.3.1) 

^_3GM/i^V_ (2) 



\ C , 

where we are assuming a circular orbit {6 = 0, and thus 9 = n) of semimajor axis a about a planet 
with mass M. Forced libration due to orbital eccentricity is considered in Section [5| 

This is simple harmonic motion, with a solution of the form 

i/j = a cos{ujot + (3) , (3) 
where a and f3 are constants to be determined by satisfying initial conditions and 



3GMB-A ^(B-A) 

Given our assumption of an instantaneous orbit swap (see Section |2]), the initial conditions are 
'i/'(t = 0) = ni — n2 and 'i/'(t = 0) = 0, resulting in 

ip = smcJoi- (5) 

Wo 
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Using orbital parameters from Jacobson et al. (20081 and shape parameters from Section |6] listed 
in Table |2] we find the amplitude of the free libration in this case is 0.022° for Janus and 0.046° 
for Epimetheus. 



4. Free Libration: The Damped Case 

Dissipation within the interior of the moon creates a frictional damping torque on the libration. 
Parameterizing the friction as b, we now have the damped harmonic oscillator: 

'ip + bip + wgV = 0. (6) 

The only case in which libration proceeds is the under-damped case {b/2uJo < 1), for which the 
general solution is 

^ ^ g-w/2 (cKCoswt + /3sina;f) , (7) 
where a and f3 again are constants to be determined, and 

Applying our initial conditions, as before, in the case of an instantaneous orbit swap, we obtain 




^jj='-^^^e-''/'smo,t. (9) 



What is the relation between b and the tidal dissipation parameter Ql Following the derivation 



of Murray and Dermott (1999 p. 161), the peak energy stored in the oscillator per unit mass is 

Eo= f u^l^di, = \ujIj^^ (10) 
JO ^ 

for instantaneous amplitude is/, and the energy dissipated per unit mass per cycle is 

AS= -6(i^w)2. — , (11) 
2 00 



giving 

n = _ 

A^; bu) 

Table 2: Parameters for free libration 



<J=^ = ^. (12) 





ni (7dy) 


na (Vdy) 


ni - 712 (°/dy) 


{B-A)/C a;o (°/dy) 


Amplitude (°) 


Janus 


518.238 


518.346 


0.108 


0.100 284 


0.022 


Epimetheus 


518.486 


518.098 


-0.388 


0.296 488 


0.046 



Note: ni and n2 are taken from | Jacobson et al. ||2008^ corresponding to intervals before and after the 



2006 orbit swap, while {B — A)/C comes from Table [6] below. 
Amplitude is i/iq, the value of i/i when t = in the undamped case (Eq.jsjl. 
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Substituting for lo from Eq. [s] we solve the resulting quadratic in (ft/wo)^- Rejecting the positive 
branch of the quadratic solution because of its unphysical behavior at large Q, we obtain 



where the latter approximation is valid for \. In practice, h k, uj^jQ holds for Q^2. 

We can now write the damping e-folding timescale, r = 2/6, as 0.40Q days for Janus, and 
0.24(5 days for Epimetheus. We again use Q ~ 100 for a solid icy body and Q ~ 10 for a fractured 
rubble-pile, but the difference in damping timescale between the two cases is now only one order of 
magnitude instead of two (as it was for the naive treatment in Section [2]), as the rigidity no longer 
enters into our equations. 

The initial libration amplitudes in the damped case are equal to the values listed in Table [2] 
divided by the quantity y^l — IP /Auj^^ which is very close to unity for the values of Q we expect 



(see Eq. 13 ) 



It is worth considering whether the energy dissipated in the damping of the free libration 



might be a significant heat source for the moons. From Eq. 10, using values for ujq and given in 
Table [2] and using the mean radius R (Table [l| to convert the angular amplitude into one in units 
of distance, the initial energy per unit mass stored in the oscillator is ~ 10~^ J/kg for Epimetheus, 
and ten times less for Janus. Given the idealization of our model, as stated before, this is an upper 
limit. Even if this energy were dissipated very quickly, given the specific heat capacity of water ice, 
the moons would be heated by only ~ 10~^ K. 

A possible refinement to this model would be a more sophisticated treatment of the mean- 
motion evolution during an orbit swap, perhaps in a numerical integration or by use of a poly- 
nomial fit to the mean motion from a numerical integration. However, the current model of an 
instantaneous swap, when compared to any more realistic model, will maximize the resulting free- 
libration amplitude. Thus our calculated amplitudes should be interpreted as upper limits, which 
are sufficient given our conclusion that the free librations induced on the two co-orbitals by the 
orbit swap are not currently observable. 

The result of this calculation is that the amplitude of the free libration excited by the orbit 
swap of Janus and Epimetheus is always less than 0.1°, regardless of our choice of Q, even in the 
idealized case treated here. Furthermore, the damping timescale is only a few weeks, even in the 
case of very weak damping {Q ~ 100). Current Cassini Imaging data have neither the spatial 
resolution nor the temporal sampling necessary to detect the free libration. 



5. Forced Libration 



Rotational libration can also be forced by any orbital eccentricity of the moons. We still 
define ip as the angle between the moon's long axis and the direction towards Saturn, and we 
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moon ^ 

= 2e*sin(nt) 


\ 


/ 









a>o > n 

/^f\^\'^ " 2e*sin(nt) | 




/ 









Fig. 3. — The orientation of the moon's long axis is measured from the direction towards Saturn (-0) 
and from the direction towards the empty focus of the moon's orbit (7). The latter is the deviation 
from synchronous rotation for e <C 1. All angles are positive in the counter-clockwise direction, 
so both and 7 are negative in the upper panel. The optical libration angle (\) is determined by 
the eccentricity and orbital phase. The natural frequency wq, which is determined by the shape 
parameter [B — A) jC (Eq.|4|, determines whether the moon's long axis points away from or past the 
planet, as seen from the empty focus of the moon's orbit (towards which a synchronously-rotating 
moon would face). The transition at uq ~ ''^ occurs when {B — A)/C ~ 1/3. 



now additionally define 7 as the angle between the moon's long axis and the empty focus of the 
moon's orbit (Fig. [s]). The former is more important when considering tidal effects, while the 
latter is the deviation from synchronous rotation (the "physical libration"), due to the well-known 
result that, to first order in eccentricity, a synchronously-rotating moon always keeps the same face 



towards the empty focus (e.g., Murray and Dermott 1999, p. 44). The optical libration angle is 
(f) = 2e sin nt = ^ — ip. 



This problem has the solution (e.g., Murray and Dermott 1999 Eq. 5.123) 



2e 



1 - {n/uoY 



sin nt. 



(14) 
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which is equivalent to 

2e , , 

^=7 i~v2 7Sinnt, (15) 

where e is the eccentricity and t is the time since the most recent periapse. 

If the orbit of the moon is given, the amplitude of the forced libration is then determined by 
the natural frequency wq, which arises directly from the moon's shape. Rearranging Eq. [4] we have 

:)^-(^)- 

Most moons are close enough to spherical that *C with the result that 7 approaches zero 
and approaches —4> (Fig. [3^). That is, the moon responds so sluggishly to Saturn's torques that 
it remains close to synchronous rotation, with the same face always towards the empty focus. On 
the other hand, a moon might be so elongated that ujq 3> n, with the result that approaches 
zero and 7 approaches cj) (Fig. [sja), which is to say that the moon responds so rapidly to Saturn's 
torques that it keeps the same face towards Saturn despite its eccentric orbit. 

The range < 7 < 0, which is equivalent to —cj) < -0 < 0, is forbidden. That is, the moon's 
long axis can never point between the direction towards Saturn and the direction towards the empty 



focus. We note that, due to an apparent drafting error. Figure 5.16 of Murray and Dermott (19991 
illustrates this forbidden configuration. 

The most interesting behavior occurs near the transition between these two regimes, where 



loq ~ n, which is to say that the moon's shape is such that {B — A)/C is near one-third (Eq. 16). 
In this case, small divisors cause 7 and ip to become both very negative {ujq < n) or both very 
positive {ojQ > n), as seen in Fig. |4| An extreme case of this is Hyperion, whose shape brings it 
close enough to (Jq ~ f^ so that, in combination with its relatively high eccentricity, our assumption 
of low libration amplitude is violated — as can be seen by plugging Hyperion's values, {B — A)/C ~ 



0.26 and e ~ 0.1, into Eqs. 14 and 16 Instead, Hyperion has a three-dimensional rotation state 



characterized by chaotic tumbling ( Wisdom et al.|1984 Wisdom|l987 ). As we will see, Epimetheus 



is close enough to this regime that its libration amplitude varies sensitively with its shape (Fig. |4]), 
allowing us to determine its gravitational moments with greater precision by measuring the libration 
than is achieved with direct observation of the shape. 

Preliminary calculated values of the forced libration amplitude are listed in Table [3] where 
[B — A)/C is taken from the shape parameters presented in Section pi and mean values of n and 



Table 3: Parameters for shape-derived prediction of forced libration 





n (°/dy) 


e 


{B-A)/C 


■^0 (°/dy) 


Amplitude (°) 


Janus 


518.292 


0.0068 


0.100 ±0.012 


284 ± 17 


-0.33 ±0.06 


Epimetheus 


518.292 


0.0098 


296+° °^'' 


4881^5 


_8q-i0-4 



Note: n and e are taken from |Jacobson et al.| l |2008| , while (B — A)/C comes from Table |6] below. 
Amplitude is 7^/21 the value of 7 when nt = -k 12 (Eq. |14l|. 
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Fig. 4. — Predicted libration amplitudes for 7 and ■0) calculated from Eqs. 14 and 15 with nt = 7r/2 



for the case of Epimetheus, are plotted against the changing inertial moments (shape) of the moon. 
Large amplitudes occur in the vicinity of the resonant value of the shape-derived gravity parameter, 
(B — A)/C ~ 1/3, at which the natural frequency uiq approaches the orbital mean motion (the 
forcing frequency) n = 518.292°/day. 



e are used since the variations are small. For Janus, the forced libration is larger than the free 
libration due to the orbit swap, but is still quite small at 0.33°. For Epimetheus, however, the 
amplitude of the forced libration is well into the detectable range. The calculated value is 8.9°, 
though with a fairly large uncertainty since {B — A)/C is not far from one-third (compare the 
values in Table ^ to Fig. HI) . This result will be further refined by our analysis in Section ^ 



6. Data 



Janus and Epimetheus have been imaged on a few occasions at pixel scales better than 1 km. 



and dozens of times at lower resolution, by the Imaging Science Subsystem (Porco et al. 2004 1 



on board the Cassini spacecraft (Fig. [5]). Details for images used in this paper are given in 
Tables |4] and [5j From these data, a numerical model of the size and shape of each moon has been 
constructed using established techniques ( Thomas||1993 Simonelli et al.|[l993 Thomas et al.|[T998 |. 
Stereogrammetric solution of control points is the first step, and provides relative three-dimensional 
reference points over most of the surface (for both satellites, the solutions were closed around all 
longitudes). The libration solution is part of this "boot-strap" process that is necessary for reducing 
errors in the control points and deriving an accurate shape; the libration's amplitude and phase 
are involved because they determine the satellite's precise orientation as seen in any image. We 
obtained 66 control points on janus, and 49 on Epimetheus. Limb and terminator positions are then 



Fig. 5. — Janus {left) and Epimetheus {right), from Cassini images N1537923147 and N1575363491. 
Note the different scale bars. 



used to constrain the model further, and are crucial for filling in the model where there is sparse 
coverage of control points. The accuracy of the shape models is restricted by the available views 
and resolutions, and by the limited seasonal illumination that so far leaves small areas near the 
north poles poorly observed. The estimated accuracy of a shape is calculated from the resolutions 
and orientations of the applied images. Shape models for Janus and Epimetheus were given by 



Porco et al. (20071, but are here revised with a more accurate accounting for the rotation state, as 



well as using additional recent data. The improved rotational models and additional image data, 
for the first time, allow shape models to be useful for moment-of-inertia calculations. 

Volume and mean radius are computed directly from the numerical shape models, and gravita- 
tional moments of inertia are calculated from the models under an assumption of constant density. 
An ellipsoidal fit to the numerical shape model is done only to provide the reader with a first-order 
estimate of the principal dimensions, and should not be used for moment calculations other than 
crude approximations. 

In order to compile a large suite of images, taken at different times from distinct geometries, 
into a single shape model, one must make an assumption about the rotation state. Different assumed 
rotation states can be tested sequentially, and the best one can be determined by minimizing the 
residuals (measured in rms pixels) between predicted and actual locations of the control points 
in the images. While it is often sufficient to assume a constant rotation rate (which may well be 
synchronous rotation), perhaps with the superposition of a simple libration, the rotation states of 
Janus and Epimetheus cannot be described so simply. The periodic orbit swaps change the orbital 
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mean motion, and thus not only the synchronous rotation rate but also the frequency of the forced 
libration. A free libration would further complicate the model, but we neglect free librations due 
to the conclusions reached in Section IH 

We addressed the problem by compiling a sequence of binary-PCK kernels in the SPICE 



Table 4: Observing information for images used in this paper (Janus) 



# of Sub-S/C Sub-S/C Sub-Solar Sub-Solar Phase Range 



Image Identifier 


Images 


Date 


Latitude 


Longitude 


Latitude 


Longitude 


Angle 


(km) 


1487417508 




2005-49 


0.2 


154.56 


-22.39 


52.97 


100.783 


911572 


1495308742 




2005-140 


-27.21 


219.96 


-21.81 


221.74 


5.637 


357370 


1501558704 - 58826 


4 


2005-213 


-20.44 


296.1 


-20.76 


267.0 


27.24 


923600 


1501711325 - 11501 


3 


2005-214 


11.77 


248.4 


-20.73 


102.1 


146.56 


539000 


1507877360 - 78655 


19 


2005-286 


-0.073 


114 


-20.13 


8 


105.02 


879000 


1521539019 - 39613 


3 


2006-79 


-0.29 


278 


-18.16 


230. 


50.8 


493000 


1521670216 - 70692 


6 


2006-80 


0.35 


138 


-18.15 


296 


152.26 


724500 


1524901386 - 01922 


4 


2006-118 


-0.36 


303 


-17.85 


240. 


63.9 


703000 


1524964907 - 66330 


9 


2006-119 


-0.34 


235 


-17.84 


264 


33.6 


223000 


1537919879 - 23147 


22 


2006-268/269 


21 


176 


-15.78 


227 


61.7 


149000 


1563765026 




2007-203 


0.07 


106.34 


-11.32 


99.01 


13.519 


892568 


1582238183 - 40507 


10 


2008-51 


-72.1 


95 


-8.31 


44 


70.9 


175000 


1589742955 - 43095 


2 


2008-138 


42.5 


320.7 


-6.91 


59.4 


101.1 


350200 


1590458717 - 59027 


4 


2008-147 


-73.85 


304.6 


-6.73 


33.6 


83.26 


185000 


1590461157 




2008-147 


-73.05 


317.88 


-6.73 


47.42 


83.43 


213017 


1593508083 - 09823 


11 


2008-182 


-33 - 21 


107 


-6.12 


332 


130. 


36000 


1594708747 - 09210 


2 


2008-196 


71.7 


323 


-5.94 


330 


77.8 


259000 


1602109316 




2008-281 


32.79 


85.8 


-4.86 


84.57 


37.668 


1012129 



Latitudes, longitudes, and phase angle are in degrees. 

For lines referring to multiple images, variation in each parameter is in the last significant figure. 



Table 5: Observing information for images used in this paper (Epimetheus) 





# of 




Sub-S/C 


Sub-S/C 


Sub-Solar 


Sub-Solar 


Phase 


Range 


Image Identifier 


Images 


Date 


Latitude 


Longitude 


Latitude 


Longitude 


Angle 


(km) 


1487432013 




2005-49 


-0.09 


160.66 


-22.35 


61.37 


98.552 


989916 


1490836693 - 36932 


8 


2005-89 


-0.55 


204.4 


-21.84 


322.1 


115.3 


74600 


1493544135 




2005-120 


-17.89 


43.99 


-21.85 


8.61 


33.437 


1628222 


1495309992 




2005-140 


-27.56 


139.33 


-21.88 


167.57 


26.211 


344611 


1500071960 - 72512 


6 


2005-195 


34.5 


223.8 


-21.2 


305 


95.3 


88000 


1501559484 - 59606 


4 


2005-213 


-22.77 


253.9 


-20.8 


224.1 


27.71 


833300 


1501711887 




2005-214 


13.26 


207.29 


-20.76 


69.79 


138.793 


470743 


1514191733 




2005-359 


0.65 


270.8 


-19.45 


69.55 


151.95 


452136 


1516413432 




2006-20 


0.01 


236.46 


-18.84 


81.08 


149.366 


1527858 


1521538726 - 39613 


3 


2006-79 


-0.05 


268 


-17.94 


220. 


50.9 


453000 


1524904716 - 05315 


5 


2006-118 


-0.24 


304 


-17.9 


243 


62.8 


667000 


1560620123 




2007-166 


0.78 


140.49 


-12.03 


180.66 


41.888 


1693820 


1575363079 - 64219 


14 


2007-337 


-40.0 


99 


-9.53 


34 


64.5 


40000 


1596337648 - 37691 


2 


2008-215 


29.40 


167.1 


-5.52 


158.4 


35.91 


991900 


1599590509 




2008-252 


42.79 


240.56 


-5.38 


223.4 


50.625 


859395 



Latitudes, longitudes, and phase angle are in degrees. 

For lines referring to multiple images, variation in each parameter is in the last significant figure. 
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Fig. 6. — Residuals (observed — predicted image locations) for the solution of control points for 
Janus and Epimetheus, as a function of the amplitude of assumed forced libration (see Section |5]). 
We used 66 control points in the fit for Janus, and 49 for Epimetheus. The region shaded in gray, 
where < 7 < and —0 < < 0, is dynamically forbidden. The horizontal dotted line is the 
threshold used to estimate the uncertainty in the measured amplitude. 



navigation data system ( Acton|[l996[ )|*] In contrast to the more common text-PCK kernels, binary- 



PCK kernels allow for arbitrary orientation of a given body as a function of time. Our base state 
was to have the moon keep one face towards Saturn at all times (?/) = 0), whether before, during, 
or after an orbit swap. A forced libration of any amplitude can then be superposed atop that base 
state. We are thus assuming that the frequency of the forced libration responds adiabatically to 
the changing mean motion, an assumption we justify by noting that the orbit swap proceeds slowly 
(several months) compared to a single orbit (16.7 hr). 

The results of this process are shown in Fig. [6| For both moons, we found a clear best-fit 
value for the libration amplitude, with fit residuals increasing as the amplitude changes in either 
direction. The best-fit residuals are 0.372 (Janus) and 0.382 (Epimetheus) rms pixels, indicating a 
high-quality fit in both cases. 

We used two methods to estimate the uncertainty in the derived libration amplitude. Firstly, 
we simply adopted 0.45 rms pixels as a threshold indicating significantly reduced accuracy in the 
control-point solution (see the horizontal dotted line in Fig. [6]). The control points are nearly 
all crater rims; measurements are the image locations of the center of an ellipse fit to the rim. 
On objects with well-known spins and no morphological complications, solutions typically can 
reach ~ 1/3 pixels. The overall solution here is not likely to reach this common level due to 
inclusion of some low-resolution views where large less-circular craters are of necessity included. 



Available at http://iialf.jpl.nasa.gov/naif/toolkit.htnil 
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Residuals above ~ 0.4 pixels indicate some problem with input geometry or the necessity to include 
irregularly-shaped craters. Here 0.45 pixels is a safe measure of poor input geometry, that is, 
incorrect assumed libration amplitudes. Making use of a quadratic fit to the points within the 
threshold, we find 7,^/2 = —0.3° ±0.9° (Janus) and 7,^/2 = —5.9° it 1.2° (Epimetheus). As a second 
method of estimating the uncertainty, we take the best solution, and allow half the control points 
to have the best- fit error, and the other half an error of 0.5 pixels, and use the resulting average 
rms residuals as the error estimate. This method gives results nearly identical to the first method. 

A crude check can be made by considering the surface resolution of good images that affect 
the solution. For Epimetheus, this is 0.44 km/pixel, so that 0.45 pixels corresponds to 0.2° on 
the surface of a body of mean radius 58.1 km. Likewise, for Janus, the result is 0.3°. Thus the 
estimated uncertainties from Fig. |6] are reasonable and conservative. 

The parameters of the best-fit shape model are given in Table [6j The measured values given 
here for mean radius R and moment of inertia ratio {B — A)/C have been retroactively used for 
the preliminary calculations in the preceding sections (Tables [T] through |3]). 

The moment of inertia ratio {B — A)/C is taken directly from the numerical shape model. 
Though we give axis lengths (a, b, c) of an ellipsoid fit to the numerical shape model, the moment 
ratio estimated from the ellipsoid, {B — A)/C ~ (a^ — b'^)/{o? + 6^), gives a highly divergent value 
in the case of Epimetheus, whose shape in fact deviates significantly from an ellipsoid. We estimate 
uncertainties in the model moments by adding and subtracting bulges to the shape model, consistent 
with constraints from image data, for six comparison models. The maximum and minimum [B — 
A) /C are then used as estimated bounds on moment ratios. 

We are now in a position to compare the measured rotation state to that expected from 
calculations based on the shape model. With our measured values of the amplitude 7^/2 (that is. 



the value of 7 when nt = '/r/2; see Eq. 14) and the known orbital eccentricity and mean motion, we 
re-derive the moment of inertia ratio. 

For Janus, the measured libration amplitude implies {B — A)/C = 0.09l^Q Qg, which is consistent 
with, but clearly less precise than, the shape-derived value. However, it is interesting that, despite 
the large uncertainty, the measured 7^/2 = —0.3° is very close to the shape-derived prediction 
of —0.33° (Table [3|. This would tend to validate the assumption of constant density, which was 
employed in order to derive gravitational moments from the shape model. On the other hand, we 
note that Janus' sub-Saturn point (when nt = 0) is offset from the axis of the minimum moment 



Table 6: Parameters for best-fit numerical shape models 







B/MR? 


C/MR^ 


{B-A)/C 


a (km) 


b (km) 


c (km) 


R (km) 


Janus 


0.360 


0.407 


0.470 


0.100 ±0.012 


101.5 ± 1.9 


92.5 ± 1.2 


76.3 ± 1.2 


89.5 ± 1.4 


Epimetheus 


0.328 


0.469 


0.476 


296+° °^^ 


64.9 ± 2.0 


57.0 ± 3.7 


53.1 ± 0.7 


58.1 ± 1.8 



Principal moments of inertia (A, B, and C) are calculated under an assumption of constant density. 

Mean radius R, the radius of a sphere of equivalent volume, is derived directly from the numerical shape model. 

Axis lengths a, b, and c come from an ellipsoidal fit. 
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of inertia (A), by 5.2° it 1° in longitude and 2.3° ± 1° in latitude. This may indicate that Janus 
does in fact have some internal density asymmetries, though the latitude offset at least may be due 
simply to vertical librations (equivalent to a non-zero obliquity and precession), which have not 
been taken into account in the current model. Unlike Epimetheus, Janus has a significant value of 
(C — B)/A, making it more responsive to vertical perturbations. 

For Epimetheus, the measured libration amplitude implies (B — A)/C = 0.2801qq5^, signifi- 
cantly improving upon the precision of the shape-derived value. The libration-measured value is 
in agreement with the shape-derived estimate, which has a fairly large uncertainty, so we cannot 
constrain the possibility of internal density asymmetries without further improvements to the shape 
model. Furthermore, in contrast to Janus, Epimetheus' sub-Saturn point for nt = deviates from 
the minimum moment of inertia axis by less than 1° in both latitude and longitude. 

7. SummEiry and Conclusions 

Although the free libration, the "ringing" of the moons' rotation states due to their periodic 
orbit swap, is a tantalizing target, our calculations show that it is not detectable by Cassini Imaging. 
Not only is the amplitude too small (< 0.1°), even with some idealized assumptions, but the decay 
rate is too fast (a few weeks, at most). On the other hand, our calculations indicate that the 
forced libration, which arises on any non-spherical moon with an eccentric orbit, should be barely 
detectable for Janus and quite large for Epimetheus. 

Our analysis of Cassini Imaging data, combined with a model of the forced libration, yields 

improved shape models and (for Epimetheus) the third clear detection of forced libration in the solar 
system (after Earth's Moon and Mars' Phobos). Furthermore, Epimetheus' measured amplitude 
of forced libration allows significant improvements in the precision of the shape-derived estimate of 
the ratio {B — A)/C of principal moments of inertia. 

Although the measured forced-libration amplitude for Janus is indistinguishable from zero 



Table 7: Summary of results 





Janus 


Epimetheus 


Free libration, shape-derived prediction 






Amplitude 


0.022° 


0.046° 


Damping time, solid (Q ~ 100) 


40 days 


24 days 


Damping time, fractured (Q ~ 10) 


4 days 


2.4 days 


Optical libration, 4>Tr/2 = 2e 


0.8° 


1.1° 


Forced libration, 7.^^/2 






Amplitude, shape-derived prediction 


-0.33° ± 0.06° 


QO -10.4° 
-8.9 _^42o 


Amplitude, measured 


-0.3° ± 0.9° 


-5.9° ± 1.2° 


Moment of inertia ratio {B — A)/C 






Shape-derived (assuming constant density) 


0.100 ±0.012 


0.2961°;°^? 


Inferred from measured libration 


U.U9_Q 09 


0.2801°;°?? 
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(synchronous rotation) due to its uncertainty, it is consistent with the shape-derived estimate of 
{B — A)/ C. Moreover, we have detected a permanent offset of several degrees between the shape 
model for Janus and the observed equilibrium rotation state. 

We have particular interest in comparing the orientations and/or magnitudes of the principal 
moments of inertia inferred from analysis of their rotational librations with those inferred from 
analysis of their shapes. Since the latter involves an assumption of homogeneous density, we 
account for any discrepancies (such as the offset in orientation detected in Janus) by suggesting 
that internal density asymmetries exist within the moon. These may take the form of internal voids 



(Goldreich and Toomre 


19691 


(Porco et al. 


2007 


)• 



The principal step that might be taken to improve this method would be to consider vertical 
librations (i.e., non-zero obliquity and precession). Neglecting vertical perturbations is certainly 
justified in the case of Epimetheus, which is quite prolate in shape (^ < i? ~ C), though including 
them may improve the results for Janus. 



According to a companion paper by Morrison et al. (20091, sets of parallel grooves observed 



on the surfaces of Epimetheus and Phobos, but not on moons for which significant libration can 
be ruled out, also occur on Pandora and cannot be ruled on on Prometheus. We plan to apply our 
method to these and other moons in the Saturn system, in the hopes of detecting forced librations 
that can be used to probe the interiors of orbiting objects. 
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